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Abstract 

A new density functional theory (DFT) exchange-correlation functional for the exploration of 
reaction mechanisms is proposed. This new functional, denoted BMK (Boese-Martin for Kinetics), 
has an accuracy in the 2 kcal/mol range for transition state barriers but, unlike previous attempts 
at such a functional, this improved accuracy does not come at the expense of equilibrium properties. 
This makes it a general-purpose functional whose domain of applicability has been extended to 
transition states, rather than a specialized functional for kinetics. The improvement in BMK rests 
on the inclusion of the kinetic energy density together with a large value of the exact exchange 
mixing coefficient. For this functional, the kinetic energy density appears to correct 'back' the 
excess exact exchange mixing for ground-state properties, possibly simulating variable exchange. 
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I. INTRODUCTION 



In the last decade, density functional theory (DFT) has probably become the number 
one ap_p_lied computational chemistry tool. Most exchange-correlation functionals currently 
— including the exceedingly popular Becke 3-parameter-Lee- Yang-Parr 




usedll 



(B3LYP) hybrid functional [i 



have been developed more than ten years ago. Yet these 



first-generation functionals exhibit a number of chemically important weaknesses. 



s have been proposed as successors in recent 
"]), mPW91 (modified Perdew-Wang 19910), 
VSXC (Van Voorhis-Scuseria exchange-correlation 1(3]). PBEO (PBE with 25% exact ex- 



A number of second-generation functiona 
years, such as PBE (Perdew-Burke-Ernzerhof 



change added 11]), PKZB (Prr(lo\v-Kurtli-Zu|)au-Blalia.| 12] ). 097 (Becke 1097 13]). B97-1 



(Becke 1997 variant 1 
Tozer-Handy) family 



1 1 }J39£2 (Becke 1997 variant 2[19|), the HCTH (Hamprecht-Cohen 



[ llQ, OPTX/OPTC (optimized exchange/correlation 3, Q ) , 
r-HCTH and its hybrid jsol] . TPSS ( Tao- Per dew- S t aroverov- S cuseria ^\ ) . and B98 (Becke 



1998 |22|). Two important avenues that have been pursued in this direction are: (a) inclu- 
sion of new variables explicitly dependent on the Kohn-Sham orbitals, such as the kinetic 
energy density r; (b) treatment of DFT as a parametrized 'semi-ab initio' method 2^. 

In particular, the most commonly used hybrid functionals all involve parametrization 
against experimental data to some extent: For instance, the mixing parameters in B3LYP 
were fitted against the binding energies of the molecules in the G2-1 set [24]. Becke's subse- 
quent 1997 functional^! involves a power series expansion in the reduced density gradient, 
parameters of which were fitted against the binding energies, ionization potentials, and pro- 
ton affinities in the G2-1 set plus some additional data. An expansion of order m requires 
3m + 4 parameters: for m = 4 and higher the fit became unphysical, and even m = 3 
appeared to exhibit signs of over-fitting. Erring on the side of caution, Becke chose m = 2, 
involving 10 parameters. 

The issue of over-fitting may be substantially alleviated by (a) admitting high-level ab ini- 
tio data to the 'training set'; (b) considering data other than energetics, such as gradients at 
the equilibrium geometry (a convenient metric for geometry errors) and numerical exchange- 



correlation potentials extracted from high-level ab initio e 



Zhao-Morrison- Parr 25] procedure. Tozer and Handy 



26 



ectron densities by means of the 



27 



followed this route, with 

o, lly derate success, as their tactical fc™ was "poody caoseufl. Greater success 
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was achieved by applying this more elaborate refinement procedure to Becke's functional 
form, yielding the orig inal HCTH/93 (Hamprecht-Cohen-Tozer-HandyQ, Q) 

functional. 

Besides this m = 4 'pure DFT' functional, a reparametrization of Becke's hybrid functional 
(denoted B97-1) was obtained by-product. 

The general performance of HCTH/93 was markedly superior to first-generation GGA 
(generalized gradient approximation) functionals such as BLYP (Becke-Lee- Y ang -Parr Q| ) . 



311 ]: in general, 



However, it did not measure up to expectations for weakly bound systems [30, 
the absence of anions in the training set meant that the functional was not adequately 
parametrized for the outer regions of the electron density. In an attempt to remedy this 
problem and generally increase the versatility of the functional, it was reparametrized (la. 
to ever larger sets of data including anions and weakly bound systems, culminating in the 
HCTH/407 functional 3- 

The numerals in the notation refer to the number of chemical 
systems used in the training set. 

Attempts to include additional variables such as the kinetic energy density r into the 
functional have been reviewed in the introduction to Ref. This latter paper represents 
an attempt to enhance HCTH with r-dependent terms. One intriguing question here is the 
actual role played by the kinetic energy density. Becke has argued [32J that local density func- 
tionals that include the kinetic energy density r can simulate delocalized exchange and thus, 
inclusion of r will improve generalized gradient functionals. Van Voorhis and ScuseriafloT] 
presented evidence that such a form might be capable of simulating exact exchange. Boese 
and Handyj^jJ have similar experience, but conjectured that r does not merely simulate 



exact exchange, but rather variable exact exchange [^J . 

This brings us to the subject of this paper. One well-known deficiency (see e.g. (33), 
) of most existing exchange-correlation functionals is their poor performance for reaction 
barrier heights: These are often significantly underestimated, and it is indeed not unusual 
for DFT not to find a barrier at all. Durant|33l noted that BH&HLYP (Becke half-and- 

), a hybrid functional with 50% exact exchange, 
was the exception to the rule in that it predicted transition state geometries and barriers 
reasonably well, despite being inferior to B3LYP (20% exact exchange) for thermochemical 
properties. Later, Truhlar and coworkers [36] and Kang and Musgrave|37l] both proposed 
reparametrizations of existing functionals (modified Perdew-Wang 9] in the former, B3LYP 
without Becke exchange in the latter case) for barrier heights: both functionals involve 
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a very high fraction of exact exchange, 42.7% in the case of mPWIK (modified Perdew- 



Wang for kinetics 



36(), 55.7% in the case of KMLYP (Kang-Musgrave-Lee-Yang-Parr[37(). 



For both functionals, the improved performance for reaction barrier heights comes at the 
expense of seriously degraded performance for ground-state properties: in particular, errors 
in atomization energies and geometries exceed those for more conventional hybrid functionals 
with exact exchange in the 15-25% range by a factor of two to three . In the case of 
KMLYP, an empirical G2-style 'high level correction' was proposed, but this is strictly a 
stop-gap solution as it will not affect properties other than equilibrium energetics. 

For researchers involved in the computational study of reaction mechanisms (such as our 

fin 

own group, see e.g. [39, 40]) the present situation is very frustrating — especially when a 
Pletta a rf oo mp e«- St ep pathway is t o be co— Q - as no single functional 
can be relied upon to get both relative stabilities and barrier heights right. 

One possible route to such a functional might be variable exact exchange. If the kinetic 
energy density were truly capable of simulating variable exact exchange, however, then a 
fixed-exact exchange functional with r-dependent terms might afford a computationally 
much simpler solution. In the present paper, we will show that this is indeed the case, and 
shall propose a new general-purpose exchange-correlation functional to be known as BMK 
(Boese-Martin for kinetics). 

Before discussing the development of the new functional, however, we would like to men- 
tion a couple of alternative approaches to tackling the aforementioned problem of transition 
states. These are connected to attempts by various groups at creating 'third-generation 
functionals'. None of these functionals are currently accurate enough to replace any of the 
second-generation functionals. 

Self-interaction-free or -corrected functionals constitute the first possibility. Here, an ad- 
ditional orbital-dependent correction is introduced[41], causing the energy to lose invariance 
to rotations of occupied orbitals. Solving these equations self-consistently is both time- 
consuming and difficult [^sj . Furthermore, no gradients or higher derivatives with respect to 
the geometry are available, which greatly reduces their usefulness in practical applications. 

The second novel approach is to use functionals with so-called local exact exchange. How- 
ever, such functionals either perform much less well when describing ground-state properties 
than the known functionals 43], or have various difficulties describing spin splittings 44]. 
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II. COMPUTATIONAL DETAILS 

In order to obtain additional reference data beyond those in the HCTH/407 data set 11611 . 

nn 

the Wl and W2 (Weizmann-1 and -2) computational thermochemistry protocols |45j,|46j] were 
applied. Basically these methods represent approximations to the CCSD(T) infinite basis 
set limit: They generate average errors for total atomization energies (TAEs) in the kJ/mol 
range, and should therefore definitely be accurate enough for our purpose. The calculations 
were carried out by means of the MOLPRO2002 electronic structure package and a driver 



noc 
4c 



script written in Perl 

The parametrization of the functional itself was carried out using a collection of ad 
programs and scripts written by one of us (ADB): modified versions of both CADPAC 
and Gaussian 03 [50] were used as DFT engines. 

Both are finite basis set DFT codes. The HCTH family of functionals was parametrized 
using a triple- zeta plus double polarization (TZ2P) basis set. We discovered early on in 
the present study that this basis set is insufficiently saturated for our purposes: we indeed 
found that in the case of the HCTH and r-HCTH functionals, some basis set superposition 
error has been absorbed into the functional. On the dependence of DFT parametrization 



on the basis set used, see Ref.|38|. After some convergence studies, we decided on the 
following combination of basis sets: Jensen aug-pc2 (augmented polarization consistent-2) 
basis setj51] for {H,C,N,0,F}, TZ3Pf+diffuse (triple-zeta plus triple polarization plus dif- 
fuse, including an additional high-exponent d function as known to be required for many 



second-row compounds 



52]) for {B, Al, Si, P, S, CI}, the Wachters-Hay basis set 53] for 
first-row transition metals augmented with 2 fig polarization functions taken from the Ap- 

n n 

pendix to Ref.|54|, the SDB-aug-cc-pVTZ basis set 54] for third- and fourth-row main group 



elements, and the 6-311+G(3df) basis set 55J] for remaining elements. We are confident that 
this combination of basis sets is sufficiently close to the Kohn-Sham basis set limit for our 
purposes. 

III. TOWARDS A NEW FUNCTIONAL 

We shall try to develop a functional with the following design goals: 
• Mean absolute error below 2 kcal/mol for reaction barrier heights; 
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• Smallest mean absolute error possible for atomization energies, geometry and harmonic 
frequencies of ground-state molecules; 

• No separation between exchange and correlation functionals, in light of error cancel- 
lation between the two terms; 

• A mathematically 'simple' functional in order to facilitate insight; 

• A power series expansion, yielding the maximum number of adjustable parameters 
(usually on the order of 10-20); 

• A functional that obeys some of the more simple scaling conditions and furthermore 
is capable of modeling the exchange holej^. Isfil]: 

• As diverse and balanced a parametrization set as possible, including transition metal 
(TM) complexes and hydrogen-bonded systems. 

We will consider two different functional forms, one the hybrid HCTHI141 form (of which 
and B97-1 141 are special cases), the other the hybrid r-HCTH form. The exchange- 
correlation energy (with a as the mixing coefficient) is defined as: 

Exc = E x ,i + E x ,n-i + E c + axE HF (1) 



E 



x.i 



E / e£ DA M9x.M)* (2) 

/-r 



9Xa,l = ^C X a,l,iU X a ( 3 ) 



i=0 



E 



X,n-l 



E / e L xl DA (p«)9x«,n-i(sl)fx«Mdr (4) 

a J 

fxa(w a ) = w a -2(w a ) 3 +(w a ) 5 (5) 

g(ftr*) 2 /34/3 1 

W - = 3, fi 2W3 ( 6 ) 

M 

9Xa,n—l = 2—i C Xa,n-l,i u Xcr (7) 
j=0 

U X a = JXaSl(l + JXaSlr 1 (8) 

Ec = E Eccrcr + Eca/3 (9) 
a 

Ecaa = I e£° A (p a )g Cc rM)dr (10) 
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M 

9Caa = Y. C C™,i U C*v (H) 
i=0 

U C aa = lCaaSl{l + ICaaSlY 1 (12) 

Ecap = J e Ca p(pa, Pp)gcaf3(s 2 avg )dr (13) 

M 

9Caj3 = 2 , C Ca0,i U Ca8 

(14) 

i=0 

= 7Ca/3sL g (l +7Ca/3S^) _1 (15) 

Here, s is closely related to the reduced gradient: 

si = |Vp,| 2 p CT - 8/3 (16) 



sl 9 = \(*l + 4) ( 17 ) 



The coefficients 7 are non-linear and fixed to: 

1Xa = 0.004 (18) 

ICaa = 0.2 (19) 

7Cq/3 = 0.006 (20) 

(21) 

and the kinetic energy density r is defined as: 

r a = E(V0i.) 2 (22) 

For the B97 form, Ex, n -i will be zero. The B97-1 functional jl^ has 12 linear coefficients with 
the sum over M in equations 3, 11, and 13 going to M=2. The r- HCTH hybrid functional^ 
had 16 linear coefficients with the sum over M in equations 3, 7, 11, and 13 going to M=3. 
In both cases, the sum was truncated when no appreciable further improvement was seen 
from the inclusion of higher-order terms. The main improvement in the r-HCTH hybrid 
functional obviously came from the higher order terms which included the kinetic energy 
density. 

In a first 'numerical experiment', we varied the amount of Hartree-Fock exchange (the 
variable a in equation 1) in each of these functionals in small increments between and 



50 %, and 'self-consistently' re-refined all other parameters to 



which is similar to the G2-2 dataset of Pople and coworkers 



the HCTH/147 dataset 



24 



57. 



The RMS errors 



for both 'meta-functionals', as a function of the percentage of Hartree-Fock exchange, are 
plotted in Figure [TJ 

We can make the following observations. First of all, both curves exhibit minima in the 
20% range, consistent with experience. Second, the improvement afforded by the kinetic 
energy density terms is not nearly as large as one might have expected. Third, the main 
difference occurs far away from the minimum: the 'meta-functional' including r deteriorates 
much less rapidly than the one without, both as exact exchange is reduced to zero and as it 
is cranked up into the 50% region. 

It may be instructive to consider the evolution of the parameters — particularly the 
zeroth-order linear coefficients of equations 3-14 which contribute most to the energy — as 
a function of 'exact exchange' admixture. All these coefficients, save the second non-local 
one, will be unity for the uniform electron gas. For the B97-like form, the zeroth-order 
exchange-coefficient cx a ,i,o varies from 1.1 (at 0% HF exchange) to about 0.4 (at 50% HF 
exchange). For the r-HCTH-like form, at 50% exact exchange mixing, this coefficient stands 
at 0.5, or 1.2 higher than without the terms including the kinetic energy density. As we would 
have expected, the cxa,n-i,o coefficient changes from near-zero (0.001) at 0% HF exchange 
to -0.38 at 50% HF exchange. Hence, the terms including the kinetic energy become much 
larger, apparently 'correcting back' the exact exchange mixing. Furthermore, the zeroth- 
order correlation coefficients are affected a lot by adding exact exchange: For example, 
ccw,o actually becomes negative. For the B97 form, these coefficients change from 0.4 to 
-2.8. The r-HCTH form, however, behaves completely differently, with cxaa,o remaining 
around 0.4. For the B97 form, this effect is cancelled by the cx a /3,o coefficient, going from 
0.5 to 1.8. Again, the r-HCTH form functionals coefficient does barely change from 0.7 
to 0.9. As a conclusion, we note that the terms depending on the kinetic energy density 
have to be responsible for the relatively shallow 'r-including' curve in Figure H Moreover, 
changes in the coefficients are by and large limited to the exchange part, while for the r-free 
(B97) form, both exchange and correlation coefficients are quite sensitive to variation in the 
amount of Hartree-Fock exchange. Thus, in short, introduction of Equations 4 and 5 into 
the functional makes it capable of correcting back for excess exact exchange. 

This immediately begs the question whether this property of the r-HCTH hybrid form 
also applies to reaction barrier heights and transition state geometries. Furthermore, what 
would be the behavior of the various terms and parameters as a function of the percentage 



8 



of Hartree-Fock exchange? 

For this purpose, we included such transition states into our 'training' or 'fitting' set. At 
this stage, we employed the set of reactions proposed by Truhlar and coworkers for which 
accurate experimentally derived barrier heights (including tunneling and other quantum 
corrections) are available. To this we added two identity Sj^2 reactions for which accurate 
barriers were previous determined by W2 theory j^]. For Truhlar's reactions, we performed 
Wl and W2 reference calculations at the QCISD/MG3 geometries given in the Supporting 
Information to Ref. 5j|. Results are summarized in Table |l] Overall, the agreement between 
theory and the values derived from experiment is quite good: any discrepancies would be 
primarily due to post-CCSD(T) correlation effects, or uncertainty in the measurements or 
quantum corrections. The transition states of the reactions H 2 + F and OH + CI exhibit 
severe nondynamical correlation, and CCSD(T) data would not be very meaningful. In the 
case of the H 2 + F reaction, we prefer to use the barrier of 2.31 kcal/mol obtained by Kallay 

case of the OH + CI reaction, we ended up employing the 
experimental value of 9.8 kcal/mol. Altogether 44 barriers (22 forward and 22 reverse) are 
reported in Table HI although obviously, a method that would get both the forward barrier 
and the reaction energy right would automatically give the right answer for the reverse 
barrier — so either the reverse barrier height or the reaction energy would be redundant 
data. For determining the functional, we included these 22 forward reaction barrier heights 
plus the reaction barriers of CI + CH 3 C1 CI + CH 3 C1 and F + CH 3 F -> F + CH 3 F. 
Since only some of the molecules in Table |J were included in the HCTH/147 set, we only 



use 14 o 
energies 



the 19 reaction energies in the small set. These were augmented with nine reaction 



62j that were used by Schiitz et al. to assess local correlation methods[63|. 



If we were to exclusively consider errors in atomization energies, it is at least in principle 
possible that just a particular type of bond is poorly described, and the error statistics 
obtained are not representative for reactions involving formation and/or breaking of other 
bond types. Adding reaction energies to the training set at least allows us to verify that the 
various molecules involved at least yield correct energetics relative to their molecular 'next 
of kin'. Furthermore, typical DFT studies are probably more likely to focus on reaction 
energies than on atomization energies. 

The original HCTH/147 setQ] was thus enlarged by 22 reaction barrier heights, 23 reac- 
tion energies, the dissociation energy of the NH 3 and (H20)(NH 3 ) dimers, and some potential 
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energy points of the NH3 dimer surface 



641 ]. The latter two properties were given low weights, 



while the reaction energies were assigned ten times the weight for the atomization energies, 
thus taking into account that said reaction energies are an order of magnitude smaller than 
typical atomization energies multiplied by the molecularity (number of molecules involved) 
of the reaction. 

Once again, all parameters were 'self-consistently' refitted for every value of the exact 
exchange coefficient, and this analysis was repeated for various relative weights assigned to 
the barrier heights in the 'penalty function'. For the r-free form, profiles of the RMS errors 
for all energetics other than barrier heights are plotted in Figure El (every curve representing a 
different weight factor for the barrier heights), while the corresponding profiles for the barrier 
heights themselves can be found in Figure EH With a zero weight for the transition states, 
the expected single minimum around 20% Hartree-Fock exchange is obtained. However, as 
soon as the barriers are even given a unit weight, a second minimum appears around 40%. 
As the weight for the transition states is increased, both minima persist at least through 
Wts = 50. Note that the curves are not perfectly smooth since the final parameters are 
somewhat dependent on the initial parametrization chosen. 

For a functional to be useful for transition states, barriers have to be predicted quite 
accurately — especially for reactions with low barriers which might otherwise be spuriously 
rendered barrierless. A 2 kcal/mol RMS target accuracy would, according to Figure EH 
require Wts — 50. Using this relative weight, we have plotted RMS error profiles for both 
barrier heights and other energetic properties as a function of the percentage of Hartree-Fock 
exchange in FigureEJ for both the r-free and the r-including forms. Parameter optimizations 
for the latter took the final parameters for the former as initial guesses. 

Note that the minimum around 20% for the non-barrier height energetics is now very 
shallow, and that both functional forms exhibit a deeper second minimum around 40%. For 
the barrier heights, there only appears to be a single minimum, at 42% to be precisely. 
This value was retained for the final functionals. 

The B97-1 like r-less form was frozen at this stage, and will henceforth be denoted B97-K 
(Becke-97 for 'kinetics'). 

The r-dependent form was then subjected to its final refinement. On the one hand, we 
expanded the 'training set' to an enlarged version of the HCTH/407 setQ, encompassing 
464 energetic properties and 4008 Cartesian gradient components. This set includes nearly 
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the full G3-99 set of molecules minus some very large organic systems, plus several dis- 
sociation energies of hydrogen-bonded complexes and stationary points on the ammonia 
dimer potential surface taken from Ref . 6^ , and several ligand dissociation energies of tran- 
sition metal complexes — in addition to the reaction barriers heights already discussed. On 
the other hand, we replaced TZ2P by the large basis set described in the Computational 
Methods section. Full details of all the data are available in the Supplementary Material 67]. 

The final functional thus obtained is denoted BMK (from the present authors' initials and 
the K of 'kinetics'). The parameters of both B97-K and BMK are given in Table ITT1 Some 
of their values are somewhat surprising if we consider the coefficients of these functionals 
when no barrier heights are included, as discussed above. We will first discuss the GGA 
(B97) functional form. With a larger fraction of exact exchange, the zero-order correlation 
coefficients ccw,o and cca/3,0 will 'try to correct back' the large amount of exchange mixing, 
together with a smaller local exchange coefficient cx a ,iocai,o- Introducing transition states into 
the fit to obtain B97-K will reverse this effect, since this large amount of exact exchange is 
needed. 

The r-HCTH form, however, behaves the opposite way. Here, with a larger fraction of 
exact exchange, the non-local exchange coefficients can now 'correct back' the exchange, 
yielding positive correlation coefficients ccw,o and cc a f3,o- 

When introducing the transition states into the set to obtain BMK, these non-local 
coefficients are now needed for something else, and the zero-order correlation and exchange 
coefficients look similar to the B97 functional with a larger fraction of exact exchange. 
While having a small effect on the other coefficients at the minimum in Figure Q (see Table 
II of reference j3| for comparison), the kinetic energy density has a very large effect as 
soon as the amount of exact exchange is not around 15-25%. This explains the previously 
reported observation that functionals employing the kinetic energy density are capable of 
simulating exact exchange. It also implies that with large values of the exact exchange 
mixing coefficient, the benefits derived from incorporating the kinetic energy density and 
similar variables into DFT functionals will be even larger. 

We shall now finally assess the new functionals for a variety of properties, in particular 
as to how they measure up compared to other functionals designed either specifically for 
barrier heights or for equilibrium energetics. 
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IV. ASSESSMENT OF DENSITY FUNCTION ALS 



As our final parametrization set of 464 systems is larger and arguably more diverse 
than any previously used, we shall use it for comparison with other functionals as well. 
We have considered (a) the GGA functionals BP86 (Becke exchange^ with Perdew 1986 
correlation Q), BLYpflQ, HCTH/407 lfj, and PBE (Perdew-Burke-Ernzerhoflsl); (b) the 

fc) the hybrid functionals B3LYP 6|, 
r-HCTH hybridU, KMLYp|3J and mPWlK[3ij. 
We employed the same large spdf basis set for the assessment as was used for the final 
parametrization cycles of BMK. Summary error statistics for all functionals, broken down 
by type of molecule or property, are presented in Table 11111 

For the non-hybrid functionals, the results of the previous investigations with the TZ2P 



meta-GGA functionals r-HCTH 3 and VSXCrJ; ( c lJ: ] 
Q, B97-lQ, B97-20, PBEo|sJl]|, r-HCTH hybridQ 



and cc-pVTZ basis setsj 



are corroborated: For all molecules, the r-HCTH functional 
yields the lowest errors, followed by the HCTH/407 and VSXC functionals. Particularly 
BP86 and PBE yield very large errors when considering geometries and atomization energies. 
Interestingly, this pattern is not necessarily carried over from the TAEs to the reaction 
energies: for instance, VSXC has half the RMS error for atomization energies of BLYP 
but both functionals perform similarly (somewhat poorly) for reaction energies. As to the 
geometry gradient errors which correlate with geometry errors, see also Refs. Q, Q, the 
best performance appears to be put in by those functionals that perform best for TAEs (i.e. 
the meta-GGAs and HCTH/407). 

HCTH/407 and both meta-GGA functionals in fact are overall competitive with the 
very popular B3LYP hybrid functional, in turn surpassed by B97-1 and B97-2. Interestingly 
(although perhaps fortuitously), all hybrid functionals perform worse than the (meta-)GGAs 
for transition metal ligation energies, although this problem is greatly mitigated in r-HCTH 
hybrid. For our dataset, the best hybrid functionals (B97-1 and B97-2) overall represent an 
improvement over (meta-)GGAs for TAEs, but not significantly so for ionization potentials 
and electron affinities. As expected, the transition-state optimized mPWIK and KMLYP 
functionals perform rather poorly for the thermochemical properties, although this would 
be mitigated somewhat in the case of KMLYP by applying their recommended 'high-level' 
additive correction. For geometric parameters, the RMS error of mPWIK is twice, and that 
of KMLYP thrice, that of other functionals. Interestingly, mPWIK performs reasonably 
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well for the reaction energies. 

Our B97-K functional fares better than mPWIK, with errors for equilibrium energetics 
generally in between mPWIK and B97-1, and errors for geometries comparable to B3LYP or 
r-HCTH hybrid (cut in half compared to mPWIK). But it is for BMK that performance is 
particularly pleasing: its performance for most energetic properties is comparable to the best 
hybrid general-purpose functionals, despite its high percentage of Hartree-Fock exchange. 
Indeed, for transition metal complexes it has the second best performance of all hybrid 
functionals, after r-HCTH hybrid. For geometries, its RMS error is somewhat deteriorated 
compared to B97-1 or to its closest kin, r-HCTH hybrid, but still much less so than that of 
mPWIK or KMLYP. The performance of BMK for TAEs of neutral molecules is especially 
pleasing. Error statistics for ionization potentials and electron affinities are comparable 
to B3LYP. Overall, BMK can hold its own with general-purpose hybrid functionals for 
equilibrium energetics and geometries. 

Interestingly, the improvements previously reported j^j] for the r-HCTH hybrid functional 
compared to B97-1 now vanish because of the larger basis set (compared to TZ2P used in 
the earlier work). This suggests again that some basis set incompleteness error may have 
been absorbed in the parametrization of the HCTH family j^. 

For reaction barrier heights, the general-purpose hybrid functionals represent an improve- 
ment over all (meta-)GGAs except perhaps r-HCTH. Yet they are clearly trumped by the 
transition state-optimized functionals, three of which (mPWIK, B97-K, and BMK) have 
errors around 2 kcal/mol. 

Histograms of the errors for the BMK and B97-K functional, as well as for B3LYP 
and mPWIK, are given in Figures S— 1 through S-4 in the Electronic Supplement to this 
paper [3]. The error distributions can be fitted rather well (correlation coefficients R ranging 
from 0.98 for B97-K to 0.993 for BMK) by a normal distribution A exp(— (x— /i) 2 /2a 2 ), which 
has been displayed in each graph together with its mean \i and standard deviation a. We 
can see that we did not merely remedy an additive bias (with /x shifting from -4.1 kcal/mol 
for mPWIK to -0.6 kcal/mol for BMK), but that the error distribution was narrowed by a 
factor of two (from <r=6.7 kcal/mol for mPWIK to 3.3 kcal/mol for BMK), and that BMK 
indeed compares favorably to B3LYP (/i=-1.4, <r=4.3 kcal/mol) in both regards. In contrast, 
B97-K only represents a minor improvement (/i=-3.9, <r=6.4 kcal/mol) over mPWIK. 

Overall, BMK turns in a very strong performance in Table IIII| with energy errors in 
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close proximity to hybrid functionals that employ a much smaller fraction of exact exchange 
(15-25%). The gradient and thus geometry errors are somewhat larger, and probably in the 
range of a good meta-GGA or GGA functional such as HCTH/407. 

Another important indication of the reliability of the functionals are frequency calcula- 
tions. These are especially important in order to find out if a functional has been 'overfitted'. 
In this instance, the gradients would become quite small and achieve a small gradient er- 
ror for the wrong reason: The gradients at any geometry would be small, even if far away 
from equilibrium. This would result in a flat potential energy surface and thus in too low 
frequencies 7j| . In order to ensure a balanced evaluation set (including organic molecules) 
and reliable data, we have mainly used CCSD(T)/cc-pVQZ, CCSD(T)/aug-cc-pVQZ, and 
CCSD(T)/AN04321 harmonic frequencies from the literature (we may assume that these 



values are quite accurate because of known error cancellation 
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between core-corre 
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at ion 



84, m 



and post-CCSD(T) correlation effects) 0, Q, H 
Using these reference values, we included the following molecules: C2H 2 , C2H4, C6H 6 , CCI2, 
CF 2 , CF 4 , CH 2 NH, CH 4 , CO, C0 2 , F 2 , FCCH, H 2 , H 2 CO, H 2 0, H 2 SiO, HON, HF, cis- 



HSiOH, trans-HSiOH, N 2 , ds-N 2 H 2 , iso-N 2 H 2 , trans-N 2 H 2 , N 2 0, NH 3 , SiH 4 , SiF 4 , S0 2 , 
SO3, and s-tetrazine. This makes for a total of 202 distinct harmonic frequencies: all 
computed values are given in the Supplementary Material^ |, while mean signed, mean 
absolute, and RMS errors over the whole set are displayed in Table IIV1 We investigated 
the performance of the BP86, BLYP, HCTH/407 and PBE GGA functionals, together with 
B3LYP, B97-1, B97-2, PBEO, KMLYP, mPWIK, B97-K and BMK. For technical reasons, 
the BMK frequencies have been calculated by numerical differentiation of gradients. For the 
harmonic frequencies investigated, GGA functionals generally underestimate, and hybrid 
functionals with large fractions of exact exchange overestimate, the coupled-cluster values. 
Among the GGA functionals, HCTH/407 clearly yields the lowest errors, underestimating 
the frequencies by an average 30 cm -1 . There is little difference between the hybrid func- 
tionals B3LYP, B97-1, B97-2 and PBEO, with an improvement of 17% by B97-1 compared 
to B3LYP. Our investigations 0] indicate that this improvement for B97-1 largely derives 
from organic molecules, which are described even better than with B3LYP. Generally, the 
geometry gradient errors in Table IIHI seem to transfer rather well into frequency errors for 
these functionals, with B97-K and BMK yielding errors comparable to a functional like 
HCTH/407. Again, KMLYP and mPWIK yield frequency errors as bad or even worse than 
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GGA functionals. 

Hydrogen bonded complexes serve as another commonly used test for density functionals. 
Ab initio calculations on such systems converge quite rapidly with the electron correlation 
treatment (thanks to the long-range character of the interaction studied) but very slowly 
with the basis setj^], both with and without counterpoise (stJ corrections. Since the use of 
extrapolation schemes almost eliminates the basis set superposition errors for atomization 
energies j3] and hydrogen-bonded complexes such as the ammonia dimer [64]. it is safe to 
assume that use of standard W2 theoryj^ will suffice. Further touchstones for hydrogen- 
bonded complexes include the lengths of the hydrogen-bonds and hydrogen-bonded shifts 
calculated using CCSD(T)/aug'-cc-pVTZ. The harmonic frequency shifts are obtained at 
the MP2/aug'-cc-pVTZ level. We have computed reference values for nine hydrogen bonded 
complexes which were also used in a previous 

study H, 

in addition to the ammonia dimer. 
All these values (Table |V}, with the exception of the MP2 frequency shifts (for which 
we would expect a deviation of about 15%), should be very accurate. In Table IVI[ the 
mean % and RMS % errors are defined by: mean = i Y,? ( ^ue(™k™^(!)) _ x 100 % 

-mi / value(complcx(i)) x ~ 



and RMS = y^Ef ( ^(Swff) ~ X ) x 100 % ' In this Table ' a11 interaction energies 
(MP2/a'pVQZ and DFT/aug-TZ) have been counterpoise-corrected. For the DFT and 
MP2 methods tested, counterpoise corrections are generally on the order of about 0.5- 
2.0 kJ/mol, which is sometimes still quite large relative to the interaction energies being 
considered. For the reference W2 values, all residual basis set superposition error (BSSE) is 
likely to have been eliminated. Thus, we would expect the mean error to change by +2-3%, 
so that the mean percentual error of MP2 in Table IVl will be approximately zero, and the 
functionals underestimating the interaction energy (like HCTH/407 and B3LYP) will become 
more accurate. All the computed values can be found in the supplementary material 67 1. 
Generally, it is believed that MP2 renders much more accurate results than density functional 
theory for hydrogen-bonded systems. Overall, the consistency of results of MP2 is only 
matched by one functional, surprisingly B97-K. The errors of this functional are low for all 
tested properties, from the dissociation energy to the frequency shift of the hydrogen bond. 
BMK, on the other hand, renders a mediocre performance, which looks quite similar to the 
B97-2 functional. Concerning the dissociation energies, HCTH/407 is the only functional 
giving an RMS % error of less than 10. Only the hydrogen bond distance itself is consistently 
overestimated by about 4%. The same is true for B97-1: For all properties investigated it 
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yields low errors, except for the geometry shift. Generally, hybrid functionals are more 
accurate than GGA functionals (with HCTH/407 yielding particularly good dissociation 
energies and BLYP hydrogen bond lengths), which are in turn much more accurate than 
the other two 'kinetics' functionals investigated. mPWIK shows the behavior we would 
expect from the HF method, i.e., underestimating the interaction energy, the hydrogen and 
frequency shifts by about 40%, and overestimating the bond length. KMLYP on the other 
hand shows some behavior normally associated with a functional fit to the uniform electron 
gas, such as LSDA: The interaction energy and the shifts are overestimated and the bond 
length underestimated. Overall, the new functionals show much improved performance for 
hydrogen bonded complexes, with B97-K, surprisingly, clearly yields the lowest error of all 
functionals tested. 

Finally, we will consider performance of BMK and other functionals for transition states 
outside our training set. Unlike for the latter, all transition states will be optimized with 
the functional and basis set used. The results are presented in detail in Table S-l in the 
Electronic Supplement [3]: summary error statistics can be found in Table rVTfl The source 
of the reference values was indicated in the second column of Table S-l. Entries 1-24 are the 
reactions from the training set. The remainder are a large subset of the reactions considered 
by Kang and Musgrave in the KMLYP paper pH], and finally the reactions from the recent 



paper of Houk and coworkers 



891 ] on pericyclic reactions like ring closings, hydrogen shifts 



and Diels-Alder reactions. For this class of reactions, Houk and coworkers found average 
errors ranging from 1.7 to 3.2 kcal/mol for several hybrid functionals. It should be noted 
that the average uncertainty of the experimental values is itself 0.7 kcal/mol, casting some 
doubt on the validity of overly detailed comparisons for these systems. In one case (the ring 
closing of 1,3-butadiene to cyclobutene) , the 0.2 kcal/mol stated uncertainty appears to be 
somewhat optimistic: our computed Wl barrier deviates by 1.8 kcal/mol changing by at 
most 0.1 kcal/mol when computed at different calculated reference geometries, despite the 
absence of significant nondynamical correlation. The other reactions we investigated were 
the hydrogen shifts of 1,3-cyclopentadiene and 1,3 pentadiene, again with the result that 
the experimental values might have larger deviations than their given uncertainties: For 
the former we found a difference of 2.2 kcal/mol (which is close to their CBS-GB3 number 
reported), compared to the stated uncertainty of 0.5 kcal/mol. For the latter, the Wl result 
is somewhat closer to experiment with a deviation of 0.7 kcal/mol, the reported uncertainty 
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being 0.5 kcal/mol. In all three cases, the experimental barriers are lower than our computed 
values. 

In addition, the experimental numbers of Kang and Musgrave (marked with footnote 
symbol a ) deviate significantly from the ones reported by Truhlar and co-workers for reactions 
19-24. This casts some doubt on the accuracy of the former: However, since we do not have 
any other data available for many reactions, we will still use this set for comparison. 

In Table S-l, we have tested mPWIK, KMLYP, B3LYP, B97-1, and B97-2 together with 
the two newly obtained density functionals and the spdf basis set combination described 
in the Computational Methods section. We will refrain from the use of small basis sets 
in order to calculate transition states, in order not to blur our comparison with additional 
error sources such as BSSE and basis set incompleteness effects on reference geometries. 
During pre-optimizations, we in fact found that in some difficult cases, we were unable to 
locate a transition state with a basis set like 6-31+G*. Overall, the functionals developed 
for transition states show a strong performance throughout for the barriers. The B97-1 and 
B3LYP functionals yield very inconsistent results: While some reaction barriers are actually 
overestimated, many transition states cannot be found at all, lie below the reactants in 
energy (albeit above the entry channel long distance complex: examples were previously 
reported for Sat 2 reactions P. loll]), or are underestimated by up to 10 kcal/mol (e.g. for the 
transition state of the CH 3 F + H — > CH 3 + HF reaction, # 63). For most reactions, B97-1 
does even worse than B3LYP, usually underestimating barriers by an additional kcal/mol. 
It is somewhat peculiar that the forward reactions of the type AH + H — > A + H2 show a 
lower error for B97-1 (reactions 9-15, 36, 38, 39, 53, 63-65, and 70) than for B3LYP, with 
the opposite effect seen for the reverse reactions of this type (reactions 6, 7, 19, 20, and 49). 
B97-2 usually predicts somewhat higher barriers than B97-1. 

We will now turn to the summary error statistics in Table I VI II Compared to Table IIHI 
- where we calculated single-point energies at the QCISD geometries — the RMS errors of 
all functionals are minimized further. The BMK functional yields the lowest mean absolute 
error of 1.1 kcal/mol, showing a similar performance to B97-K and mPWIK. KMLYP yields 
results between the functionals which have 20% (B3LYP, B97-1) exact exchange and the ones 
having around 42% (BMK, B97-K, and mPWIK). While B97-2 represents an improvement 
over B97-1 (as noted in Ref. 9^|), performance for barrier heights is still clearly inferior to 
BMK, as well as to B97-K and mPWIK. (B97-2 and KMLYP are somewhat more closely 



17 



matched.) 

For the fourth set of Houk and co-workers in Table IVII1 we reported error statistics to 
either ten (25-34) or eight (25-32) reactions because of the discrepancies of the various 
computed values among each other and with experiment for reactions 33-35, particularly 
the last one. For this set, B97-1 is the only functional that underestimates the transition 
barriers on average and shows the lowest RMS error. The RMS errors of BMK and mPWIK 
are worse by about 0.7 kcal/mol, followed by B3LYP, B97-K and KMLYP. 

Finally, the results obtained for our training set (reactions 1-24 in Table S-l) are further 
validated by comparing to Kang and Musgrave's experimental numbers and their CBS- 
APNO calculations. Here, KMLYP yields errors comparable to the other functionals de- 
veloped for transition states, probably because they used this set for determining their 
functional. Again, B3LYP and B97-1 show poor performance, B97-2 less so. As Kang and 
Musgrave have already reported, G2 also yields very large errors, overestimating the exper- 
imental numbers by an average of 3 kcal/mol. For transition states, G2 seems to be less 
accurate than the functionals tested. This should not surprise the reader, considering the 
heavy reliance on low-order perturbation theory — which often has problems with transi- 
tion states because of its hypersensitivity to nondynamical correlation effects — in the G2 
method. 

Summarizing this section, without the two new functionals one is in a dilemma about 
whether to compromise on the accuracy of equilibrium energetics (with transition-state 
biased functionals like mPWIK and KMLYP, compare Tables ITTT1 IIVI and IVII with Tables 
I VI II and S-l) or on that of barrier heights (with traditional general-purpose functionals). 
B97-K and in particular BMK allow users to 'have their cake and eat it too'. In fact, overall 
performance of BMK for all properties considered here is such that it can legitimately be 
considered a general-purpose functional, albeit with its applicability extended to transition 
states. 

V. CONCLUSIONS 

Prior to this work, all exchange-correlation functionals that consistently perform reliably 
for transition states and barrier heights (such as BH&HLYP, mPWIK and KMLYP) achieve 
this feat at the expense of seriously degraded performance for all other properties. We have 
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developed two new exchange-correlation functional suitable for transition states, the hy- 
brid GGA functional B97-K (Becke-97 for kinetics) and the hybrid meta-GGA functional 
BMK (Boese-Martin for kinetics). These functionals, especially BMK, were parametrized 
against very large and diverse training sets. Both functionals achieve performance similar 
to mPWIK for transition states. However, B97-K and especially BMK perform much more 
reliably for more general properties than those earlier specialized functionals. Overall per- 
formance of B97-K is intermediate between mP W1K and the very reliable B97-1 functional 
(except for perhaps fortuitously excellent performance for hydrogen bonds); the BMK func- 
tional actually outperforms B3LYP for most properties and is competitive with the best 
currently available functionals for equilibrium energetics. Its performance for geometries 
and vibrational frequencies is somewhat less good than B3LYP but markedly superior to all 
available specialized functionals for transition states. Unlike these latter, BMK can actually 
be considered a reliable general-purpose functional whose capabilities have been expanded 
to cover transition states. 

To computational chemists engaged in mechanistic studies, BMK offers a single functional 
capable of describing all aspects of the potential surface (energetics, kinetics, structures,...) 
Initial experiments for electrical properties (to be reported on in a future paper) yield promis- 
ing results as well, in part because self-interaction error issues are less acute for functionals 
with high percentages of Hartree-Fock exchange. 

The surprisingly good performance of BMK appears to hinge on the combination of 
a high percentage of Hartree-Fock exchange with terms dependent on the kinetic energy 
density. The latter appears to afford a 'back-correction' for excessive Hartree-Fock exchange 
in systems where it is undesirable, and in effect to simulate variable exact exchange to some 
extent. Some evidence exists that the same properties of BMK make it more suitable for 
describing transition metal complexes than other hybrid functionals. These observations 
should shed some light on the problem of developing both more accurate and more widely 
applicable functionals than currently available. 

In sum, we believe we have added a very useful tool to the computational chemist's 
toolbox. 
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Note added in revision 



A referee drew our attention to a recent paper by Zhao et al. [91[ in which a successor 
to mPWIK, denoted BB1K (Becke-Becke 1-parameter for kinetics) was proposed. This 
functional combines the Becke 1988 GGA exchange functional with Hartree-Fock exchange 
and Becke's 1995 meta-GGA correlation functionally. The percentage of Hartree-Fock 
exchange is empirically adjusted to minimize the RMS error over a small dataset of nine 
barrier heights, the optimum value turning out to be 42% like in the present BMK functional. 
Clearly, this functional was not as broadly and comprehensively parametrized as BMK or 
even B97-K. Furthermore, Becke himself}"!^ states: "The functional of Part IV [i.e., Ref.j^]] 
is problematic in very weakly bound systems... and is therefore not recommended. This 
will be corrected in future publications." In summary, we have reason to believe that the 
present BMK functional will be considerably more reliable and widely applicable than the 
one proposed in Ref. 
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Figure 1: Boese and Martin, Journal of Chemical Physics 
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Figure 3: Boese and Martin, Journal of Chemical Physics 
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Figure Captions 

Figure 1: Dependence on the exact exchange mixing coefficient of the RMS error for 
the B97 and r-HCTH forms (energetic properties of the 147 set, kcal/mol). 

Figure 2: Dependence on the exact exchange mixing coefficient, and the weight assigned 
to reaction barrier heights in the parametrization, of the RMS error for the B97 form 
(energetic properties of the 147 set, kcal/mol). 

Figure 3: RMS error (kcal/mol) for barrier heights with the B97 form as a function of the 
exact exchange admixture and the weight assigned to barrier heights in the parametrization. 

Figure 4: RMS errors (kcal/mol) of the B97 and r-HCTH forms for the extended 147 
set and for barrier heights as a function of the exact exchange admixture and the weight 
assigned to barrier heights in the parametrization. 
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Table I: Wl and W2 reaction barrier heights (kcal/mol) computed in this study, compared to 
experimental data 59] 



Reaction 


forward barrier 


reverse barrier 




W 1 


Wz 


exp. 


W 1 


Wz exp. 




36.80 




OO A 

38.4 


of 1 on 

36.80 


OO A 

38.4 




19.35 




18.4 


n 0/? 

9.86 


n 

8.0 


CoHk + NHo — > CoH=; + NH-? 


11.25 




in a 

10.4 


1 t no 

17.93 


1 T O 

17.8 


CoHr + OH — > C9EU 4 HoO 


O CT O 

3.52 




O A 

3.4 


on /i n 

20.49 


on ^7 

20.7 


CH? 4 Ho — > CH/i 4 H 


1 1 OO 

11.89 


1 1 on 

11.89 


10 1 

12.1 


1/1 00 

14.88 


1 /I no 1 r n 

14.92 15.0 


CHq 4- OH — > CHa 4 O 


9.00 




T O 

7.8 


1 A OO 

14.28 


1 O *7 

13.7 


CH/i + NH — > CH^ 4- NHo 

\ — J 11/1 | _L i X X ' V_.y 1 ± \ | _L i X X V 


01 r\ 

21.98 




22.1 


n cr 

8.95 


8.4 


CH^ 4 NHo — > CH? 4 NH4 


1 O no 

13.92 




-1 r r 

15.5 


1 c* nT 

16.97 


1 r? n 

17.9 


CH/i 4 OH — > CHq 4 HoO 


O OO 

6.22 




0.7 


1 n oc 

19.66 


on 

20.2 


Ho 4- CI — > H 4- HC1 

XXV ^ V_y 1 ' 11 ( 11VJ1 


b.97 


t nn 

7.90 


8.7 


5.62 


a no CT £? 

4.98 5.6 


Ho 4 F — > H 4 HF 

11^ 1 ' 11 111 


(a) 




1.8 




OO O 

33.2 


Ho 4- OH — > H 4 HoO 

I 1 / v/ii ^ 11 1L z vy 


5.29 


5.40 


5.7 


1 *70 

21.72 


01 r* c 00 n 

21.00 22. (J 


HoS 4 H — > SH 4 Ho 

ll^U | X X ' k_JXX 1 XX^ 


O AT 

3.45 


3.62 


3.6 


1 17 n 1 

17.01 


1 T 1 A ~t *7 A 

17.14 17.4 


Ho 4- H — > H 4 Ho 


9.60 


9.66 


9.6 


n /^n 

9.60 


n o g n /■* 

9.66 9.6 


CH?OH + H-> CHoOH 4 Ho 

\ — J X X < V. ./ X X 1 -L X ' \ — J J. J. / \ / XX | XX / 


9.64 




7.3 


1 C 'Id 

15.70 


IOO 

13.8 


CHq 4- HC1 -> CH/i 4 CI 

X x < 1 XX V_« X ' Y__. X X £\ \ V_> X 


1 37 

x .O 1 




1.8 


7 3Q 


7 8 


HC1 4 H -» H 4 HC1 


17.21 


17.25 18.0 


17.21 


17.25 18.0 


N 2 H 2 4 H -» N 2 H 4 H 2 


2.90 


2.81 


5.9 


40.48 


40.64 41.0 


PH 3 4 H — PH 2 4 H 2 


2.58 


2.78 


3.2 


24.58 


24.75 25.5 


NH 3 4 OH — > NH 2 4 H 2 


3.54 




3.2 


13.92 


13.2 


OH 4 CI -> 4 HC1 


(b) 




9.8 




9.9 


OH 4 H -> 4 H 2 


10.96 


10.94 10.1 


13.25 


13.48 13.1 



(a) CCSD(T) of doubtful value due to severe nondynamical correlation. Value of 2.31 kcal/mol 



61] 



taken from benchmark calculation by Kallay et al. 
(b) CCSD(T) of doubtful value due to severe nondynamical correlation. Experimental value used 

(see text). 
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Table II: Coefficients of the new functionals, where both new functionals have 42% exact exchange. 
In addition, the B97 and r-HCTH functionals when just fitting them to the 147 set (without 
transition states) are shown for comparison. 



Functional for 


transition states + energies 


energies 


Name 


B97-K BMK 


B97-K like BMK-type 


Exact exchange mixing 


42% 


45% 


c 1 = CXa lornl 


0.507863 0.474302 


0.454579 0.554540 


C2 ^J£(T TLOYL~~loCCll 


-0.192212 


-0.808213 


Cs = CCcrrr 


0.12355 -2.19098 


-2.51595 0.45226 


CA = CCnB fl 


1.58613 1.22334 


1.68137 0.89243 


CR = CXrr lornl 1 


1.46873 2.77701 


1.70057 0.11565 


L^K C. X rr T) rttl / fieri 1 1 

\J S\. U •) 1 LU 1 v , _L 


0.47394 


0.68619 


c 7 = C CW,1 


2.65399 23.8939 


12.8997 -9.17579 




-6.20977 -3.4631 


-8.9966 6.47624 


Co = CVfj lnral 1 


-1.51301 -11.4230 


-2.2741 -0.32078 


c 10 = CXcr,non—local,2 


-26.6188 


-16.0399 


Cll = CC CTCT ,2 


-3.20694 -44.3303 


-11.0840 27.6665 


Cl2 = CCq/3,2 


6.46106 10.0731 


12.6502 -36.9632 


Cl3 = CXa,local,3 


11.7167 


-0.6402 


Cl4 = £Xcr,non—local,3 


22.4891 


-5.7755 


Cl5 = CYo-o-,3 


22.5982 


-19.9035 


Cl6 = CXa/3,3 


-11.1974 


31.8788 
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Table III: RMS energy errors for the new functionals, compared with BP86, BLYP, HCTH/407, 
PBE, r-HCTH, VSXC, B3LYP, B97-1, B97-2, PBEO, r-HCTH hybrid, KMLYP and mPWIK. 



Class 


BP86 


BLYP HCTH PBE 
/407 


t- VSXC 
HCTH 


B3LYP 


B97-1 


B97-2 


PBEO r- KMLYP mPWIK B97-K BMK 
HCTHh 


All systems (404) 


21.2 


10.3 


7.9 


18.9 


7.1 


7.8 


7.7 


5.3 


5.3 


9.3 


5.5 


25.6* 


14.0 


10.3 


5.6 


Neutrals (219) 


21.8 


11.6 


8.6 


21.3 


7.6 


6.8 


8.1 


4.8 


1.6 


8.5 


5.1 


27.7* 


15.8 


11.5 


4.5 


Anions (58) 


15.2 


9.5 


7.8 


13.4 


6.4 


10.3 


9.1 


6.3 


6.1 


9.7 


6.3 


20.9* 


13.9 


10.8 


7.9 


Cations (88) 


15.1 


8.3 


7.0 


18.8 


6.3 


9.3 


5.7 


5.5 


5.6 


11.8 


5.6 


27.7* 


10.2 


7.0 


6.6 


TM complexes (5) 


11.1 


4.2 


4.7 


13.3 


4.1 


5.5 


17.0 


13.1 


15.2 


14.3 


7.2 


29.0* 


28.7 


25.0 


10.8 


TSes (24) 


10.1 


8.5 


7.0 


10.3 


7.9 


6.2 


5.0 


5.2 


4.1 


4.9 


6.1 


3.2 


1.8 


2.0 


2.0 


React, ener. (62) 


5.0 


6.4 


5.7 


5.3 


4.6 


6.8 


4.4 


4.2 


3.4 


4.3 


4.0 


8.2* 


4.8 


4.4 


3.7 


IPs (80) 
EAs (58) 


0.23 
0.22 


0.28 
0.18 


0.25 
0.24 


0.24 
0.19 


0.23 
0.18 


0.24 
0.18 


0.23 
0.19 


0.22 
0.16 


0.21 
0.18 


0.23 
0.20 


0.22 
0.17 


0.41* 
0.27* 


0.27 
0.26 


0.27 
0.22 


0.24 
0.19 


YjAgradient (4008) 


17.52 


19.01 


12.10 


16.84 


11.56 


10.96 


11.17 


10.72 


11.80 


12.33 


10.77 


30.71 


20.24 


11.61 


13.53 



All errors in kcal/mol, except ionization potentials (eV), electron affinities (eV), and sum 
of gradient components (a.u.) at their reference (experimental/theoretical) equilibrium 
geometry. 

(*) Exclusive of 'high level correction' from Ref. < 
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Table IV: Harmonic frequency errors (cm 1 ) for a variety of functionals. 



Pi i Tt pfinn a 1 


t yi o in prrnr ttipqti 


aUoUl U. LC CI 1 \JL 


RMC! prrnr 


1)1 ou 


^7 R 
-*J ( .u 


Q 

Ou. a 


7Q 1 


RT VP 


-01.0 


£9 9 
uz.z 


oo.y 


uptu /zL07 

O V. 1_ LI 1 L t\J t 


97 R 




^9 n 
uz.u 


PRF 






7^ 




9 7 
-Z. / 


91 ft 

Z± .U 


oU.o 


RQ7 1 


9 ^ 

~Z. <J 


1 R ^ 
±o.o 


9fi 3 


RQ7 9 


1 9 fi 


99 
zz. y 


^9 R 
oz.o 


PBEO 


12.5 


23.9 


34.8 


KMLYP 


86.3 


87.0 


104.6 


mPWIK 


57.6 


58.0 


72.6 


B97-K 


28.4 


30.2 


40.7 


BMK 


27.3 


31.7 


45.9 
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Table V: Reference values for the ten hydrogen-bonded complexes. 



Dimer 




Donor ArxH 


H-Bond Length 


Freq. Shift 




(kj/mol) 


(A) 


(A) 


(cm- 1 ) 


(C1H)(NH 3 ) 


34.9 


0.0434 


1.79 


-770 


(CO)(HF) 


7.1 


0.0018 


2.07 


-32 


(FH)(NH 3 ) 


52.1 


0.0317 


1.70 


-785 


(H 2 0) 2 


20.8 


0.0063 


1.95 


-170 


(H 2 0)(NH 3 ) 


26.8 


0.0117 


1.98 


-298 


(H 3 0)(H 2 0)+ (*) 


141.2 


0.2366 


1.19 


-2780 


(HC1) 2 


8.4 


0.0039 


2.56 


-72 


(HF) 2 


19.1 


0.0057 


1.82 


-138 


(OC)(HF) 


14.8 


0.0061 


2.08 


-172 


(NH 3 ) 2 


13.2 


0.0033 


2.30 


-134 



Dissociation energies by W2 theory; geometry shifts CCSD(T)/aug'-cc-pVQZ; frequency shifts MP2/aug'-cc-pVQZ. 
aug'-cc-pVQZ stands for cc-pVQZ on hydrogen and aug-cc-pVQZ on all other elements. 
(*) Symmetric structure with central H atom, better written as H5OJ. 
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Table VI: Relative errors of the functionals and MP2/aug'-cc-pVQZ for the hydrogen-bonded 
plexes. 





Dissoc. 


energy 


Ar(H-bond) 


Ar(H-bond) 


Freq. Shift 


Functional 


mean % 


RMS % 


mean % RMS % 


mean % RMS % 


mean % 


KMb /o 


BP86 


-12.5 


25.3 


80.9 


91.2 


-2.87 


A on 

4.20 


A n o 
40.3 


/inn 
49.9 


BLYP 


-14.6 


21.6 


42.7 


46.8 


-1.01 


1.97 


i n 
lb. 9 


on o 
30.2 


HCTH/407 


-7.5 


9.0 


-2.7 


26.3 


A OO 

4.28 


5.88 


1/11 

-14.1 


OO o 

23.3 


PBE 


8.6 


14.0 


76.9 


86.0 


o on 

-3.30 


A tr a 

4.54 


/in 1 
4U.1 


/i n o 
49.3 


B3LYP 


-8.5 


15.0 


29.0 


31.3 


-0.39 


1.32 


n i 
9.1 


1 ^ o 

lb. 3 


B97-1 


3.0 


6.0 


33 4 


38 2 


-1 01 


2 07 


11 A 

11.4 


lb. 2 


B97-2 


-19.6 


25.3 


18.9 


25.0 


1.08 


2.96 


4.4 


8.5 


PBEO 


2.3 


9.4 


46.1 


52.8 


-2.46 


3.26 


24.0 


29.3 


KMLYP 


32.2 


34.5 


51.4 


57.6 


-5.38 


5.66 


33.0 


4b.O 


mPWIK 


-38.0 


40.9 


-32.1 


36.3 


9.67 


10.55 


-32.9 


40.0 


B97-K 


-1.7 


6.7 


2.8 


6.0 


1.26 


1.40 


-9.2 


13.0 


BMK 


-17.8 


24.8 


32.4 


37.9 


-0.68 


1.50 


4.5 


25.2 


MP2 


-1.8 


8.0 


15.2 


19.6 


-0.33 


1.47 


(0.0) 


(0.0) 
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Table VII: Error statistics of the functionals for the various sets of barrier heights in Table S-l. 





Set Specifications 


Method 


Set 


Reference 


# TSs 


Error 


BMK 


B97-K 


mPWIK 


KMLYP 


B3LYP 


B97-1 


B97-2 


G2 


1,2 


Wl/W2,exp 24 


mean 


-0.8 


-0.3 


-0.7 


-2.4 


-4.1 


-3.9 


-2.6 




1,2 


Wl/W2,exp 24 


mean abs. 


1.1 


1.3 


1.3 


2.5 


4.1 


i.l 


3.0 




1,2 


Wl/W2,exp 24 


RMS 


1.5 


1.6 


1.5 


2.9 


1.3 


4.6 


3.5 




1,2 


Wl/W2,exp 24 


max(+) 


1.5 


2.0 


1.7 


0.9 


-1.7 


0.8 


1.8 




1,2 


Wl/W2,exp 24 


max(— ) 


-3.7 


-3.0 


-3.2 


-5.0 


-6.9 


-7.0 


-6.4 




4 


Wl,exp 


10(8) 


mean 


1.7(1.9) 


3.2(2.9) 


2.2(1.1) 


3.4(1.9) 


1.3(0.8) 


-l.O(-l.l) 0.9(-0.1) 




4 


Wl,exp 


10(8) 


mean abs. 


2.9(2.2) 


3.4(2.9) 


2.9(2.0) 


3.9(2.5) 


3.3(2.2) 


2.2(1.4) 


2.4(1.9) 




4 


Wl,exp 


10(8) 


RMS 


3.5(2.7) 


4.4(3.5) 


3.6(2.3) 


5.1(3.2) 


4.4(2.8) 


2.9(1.7) 


3.3(2.3) 




4 


Wl,exp 


10(8) 


max(+) 


6.8(5.1) 


9.9(6.1) 


6.9(4.5) 


10.2(6.1) 


10.6(5.3) 


5.1(0.9) 


8.4(3.5) 




4 


Wl,exp 


10(8) 


max(— ) 


-4.3(-0.7) 


-0.9(0.4) 


-2.0(-2.0) 


-1.2(-1.2) 


-4.5(-2.2) 


-6.0(-3.0) 


-3.1(-3.1) 




3 


CBS-APNO 


22 


mean 


0.0 


0.4 


0.5 


-0.8 


-4.1 


-4.0 


-2.6 


2.1 


3 


CBS-APNO 


22 


mean abs. 


1.3 


1.7 


1.5 


1.8 


4.3 


4.2 


2.9 


2.2 


3 


CBS-APNO 


22 


RMS 


1.7 


2.2 


2.0 


2.3 


4.9 


4.9 


3.5 


2.8 


3 


CBS-APNO 


22 


max(+) 


5.0 


5.4 


5.5 


4.8 


2.1 


1.5 


1.6 


8.7 


3 


CBS-APNO 


22 


max(— ) 


-2.8 


-2.8 


-2.9 


-4.6 


-9.2 


-8.5 


-6.7 


-1.1 


3 


exp 


32(31) 


mean 


0.7 


0.9 


1.3 


-0.8 


-3.3 


-3.3 


-1.5 


(2.9) 


3 


exp 


32(31) 


mean abs. 


1.7 


1.7 


1.4 


1.8 


3.7 


3.8 


3.0 


(3.3) 


3 


exp 


32(31) 


RMS 


2.2 


2.4 


1.9 


2.3 


1.6 


1.6 


3.5 


(4.2) 


3 


exp 


32(31) 


max(+) 


5.6 


5.9 


5.7 


3.1 


2.2 


2.6 


3.8 


(13.7) 


3 


exp 


32(31) 


max(— ) 


-2.4 


-2.5 


-0.8 


-2.5 


-11.6 


-7.9 


-5.8 


(-3.3) 



Set 1 consists of the reactions from Table HJ Set 2 are some anionic Sat 2 reactions taken from 
Ref. 6p| ; Set 3 was taken from Kang and Musgravej^l; the pericyclic organic reactions from 
the work of Houk and coworkers! 891 constitute Set 4. 
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E-PAPS Supplementary data for paper 



Will eventually be available online at (conjectured URL) 

http : //netserver . aip . org/cgi-bin/epaps?ID=JCPSA6- 121-303431 

A mirror copy will be kept at http://theochem.weizmann.ac.il/web/papers/BMK.html 



Note that in the present preprint, Table S-l has been renumbered Table VIII to avoid 
conflicts with the arXiv . org PDF hyperlinking technology 



100 



Gaussian fit: mean -0.60; std.dev. 3.33 kcal/mol; R=0.993 




Error [kcal/mol] 



Figure S-1: Error distribution of the BMK functional for the large set of 402 energies excluding 

the reaction energies. 
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Figure S-2: Error distribution of the B97-K functional for the large set of 402 energies excluding 

the reaction energies. 
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Figure S-3: Error distribution of the B3LYP functional for the large set of 402 energies excluding 

the reaction energies. 
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Figure S-4: Error distribution of the mPWIK functional for the large set of 402 energies 

excluding the reaction energies. 
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Table VIII: Performance of new functionals for barrier heights (kcal/mol) of additional reactions. 





Reference Method 


Functional Tested 


# 


Set 


Reaction 


Wl 


G2 


CBS-APNO/ cxp. 


BMK B97-K mPWIK 


KMLYP B3LYP 


B97-1 


B97-2 








/W2 




/CBS-QB3 


















1 


1 


C 5 H 8 - C 5 H 8 


36.8 






38.4 


37.0 


38.6 


36.1 


37.8 


35.0 


33.3 


33.9 


2 


1 


C 2 H 6 + NH -> C 2 H 5 + NH 2 


19.4 






18.4 


18.2 


16.9 


18.2 


16.6 


14.2 


12.7 


14.5 


3 


1 


C 2 H 6 + NH 2 -+ C 2 H 5 + NH 3 


11.3 






10.4 


11.3 


11.7 


12.0 


9.8 


8.9 


6.7 


8.2 


4 


1 


C 2 H 6 + OH -> C 2 H 5 + H 2 


3.5 






3.4 


2.9 


4.1 


4.7 


2.5 


-0.5 


-1.3 


0.0 


5 


1 


CH 3 + OH -> CH 4 + O 


9.0 






7.8 


6.3 


6.7 


7.6 


5.4 


4.6 


2.3 


3.1 


6 


1 


H 2 + CI -> H + HC1 


7.9 






8.7 


5.1 


5.3 


4.7 


2.9 


5.2 


3.9 


5.2 


7 


1 


H 2 + F -» H + HF 


2.3 






1.8 


-1.4 


-0.7 


0.5 


-1.9 


-4.6 


-4.6 


-3.0 


8 


1 


CH 3 + HC1 -» CH 4 + CI 


1.4 






1.8 


-1.0 


-0.8 


0.5 


-1.5 


-1.0 


-3.3 


-1.8 


9 


1 


H 2 S + H -> SH + H 2 


3.6 






3.6 


3.6 


4.8 


2.7 


0.9 


0.4 


4.0 


5.2 


10 


1 


H 2 + H -» H + H 2 


9.7 






9.6 


9.5 


10.2 


7.1 


4.9 


4.1 


8.8 


9.8 


11 


1 


CH 3 OH + H CH 2 OH + H 2 


9.6 






7.3 


8.3 


9.6 


6.7 


5.3 


2.6 


6.1 


7.0 


12 


1 


HC1 + H -> H + HC1 


17.3 






18.0 


17.2 


19.3 


16.3 


15.7 


12.7 


16.5 


18.3 


13 


1 


N 2 H 2 + H — > N 2 H + H 2 


2.8 






5.9 


2.6 


3.4 


1.7 


-0.1 


-0.1* 


2.5 


3.5 


14 


1 


PH 3 + H -> PH 2 + H 2 


2.8 






3.2 


3.8 


4.2 


2.1 


0.3 


0.2 


3.6 


4.6 


15 


1 


OH + H -> O + H 2 


11.0 






10.1 


10.3 


12.3 


9.1 


7.4 


4.1 


7.8 


8.3 


16 


1 


OH + CI -> O + HC1 








9.8 


7.9 


10.0 


9.0 


8.1 


5.3 


2.9 


3.4 


17 


2 


CH3CI + CI" -► CI" + CH3CI 


2.7 








3.6 


1.8 


3.6 


3.0 


-0.7 


-1.0 


1.6 


18 


2 


CH 3 F + F" ->■ F" + CH3F 


0.3 








1.8 


0.3 


2.0 


0.6 


-2.3 


-2.2 


0.2 


19 


1,3 


CH 3 + H 2 -► CH 4 + H 


11.9 


12.8 6 


11. 3 C 


12.1,8.7" 


10.0 


9.0 


9.8 


7.44 


8.7 


7.2 


8.4 


20 


1,3 


H 2 + OH ->■ H + H 2 


5.4 


6.8 6 


4.6 C 


5.7,3.6" 


4.0 


4.7 


4.8 


2.2 


0.9 


0.4 


1.8 


21 


1,3 


CH 4 + NH 2 ->■ CH3 + NH 3 


13.9 


15.0 6 


12.9 C 


15.5,12.6" 


13.7 


14.1 


14.4 


12.3 


11.4 


9.3 


10.8 


22 


1,3 


CH 4 + OH -» CH 3 + H 2 


6.2 


7A b 


4.9 C 


6.7,6.6" 


5.6 


6.9 


7.3 


5.2 


2.2 


1.1 


2.5 


23 


1,3 


NH 3 + OH — > NH 2 + H 2 


3.5 


5.1 b 


2.5 C 


3.2,1.9" 


2.7 


4.1 


5.0 


3.5 


-2.5 


-3.5 


-1.7 


24 


1,3 


CH 4 + NH -> CH 3 + NH 2 


22.0 


23.3 b 


21. 3 C 


22.7,20.5" 


21.2 


19.8 


21.2 


19.7 


17.3 


15.9 


17.8 


25 


4 


cyclobutene — > butadiene 


35.3 d 




33.7° 


33.6" 


40.4 


37.4 


39.8 


41.4 


33.0 


34.6 


35.8 


26 


4 


cis-l,3,5-hexatriene— > 




























cyclohexadiene 






28.8° 


30.2" 


31.9 


33.2 


31.3 


32.7 


31.2 


29.4 


30.1 


27 


4 


ortho-xylylene — > 




























benzocyclobutane 






25.9° 


29.8" 


29.2 


30.2 


28.5 


30.8 


28.3 


26.4 


26.7 


28 


4 


1,3-pentadiene [1,5]- 




























sigmatropic H-shift 


39. 5 e 




38.9° 


38.8" 


40.3 


42.1 


39.7 


41.2 


38.9 


36.9 


37.6 


29 


4 


1,3-cyclopcntadicne [1,5]- 




























sigmatropic H shift 


28. 2 e 




28.1° 


26.0" 


27.5 


28.8 


26.2 


27.0 


27.2 


26.4 


25.3 


30 


4 


1,5-hexadiene [1,5]- 




























sigmatropic H-shift 






33.3° 


34.8" 


36.7 


39.1 


37.9 


40.2 


26.4 


33.6 


35.3 


31 


4 


ethylene + butadiene — > 




























cyclohcxcnc 






20.7° 


21.1" 


25.4 


27.2 


23.8 


21.9 


26.4 


22.0 


24.6 


32 


4 


ethylene + cyclopentadiene — > 




























norboncne 






15.1° 


19.4" 


21.8 


23.8 


20.3 


18.2 


24.4 


19.4 


22.0 


33 


4 


2 cyclopentadiene — > 




























endo-dicyclopentadiene 






9.8" 


13.3" 


20.l' 


23.2' 


19.7' 


21.7' 


23.9' 


18.4' 


21.7 


34' 4 


cis-triscyclopropacyclohcxane— > 




























cyclonona-l,4,7-triene 






23.5° 


26.4" 


22.l' 


25.5' 


33.3' 


36.6' 


21.9' 


20.4' 


27.4' 


35 




cis-triscyclobutacyclohcxanc^ 




























cyclodeca-l,5,9-triene 








50.0" 


67.7' 


60.3' 


65.5' 


70.3' 


48.8' 


52.8' 


55.1 



(continued on next page) 



Table VIII: (continued) 





Reference Method 


Functional Tested 


# 


Set 


Reaction 


Wl 


G2 


CBS-APNO/ cxp. 


BMK B97-K 


mPWIK 


KMLYP 


B3LYP 


B97-1 


B97-2 








/W2 




/CBS-QB3 


















36 


3 


C 2 H 6 + H -> C 2 H 5 + H 2 




13.4 6 


9.4 C 


8.8" 


11.4 


13.5 


10.5 


9.0 


6.6 


10.2 


11.0 


37 


3 


CH 3 + CH 4 -> CH 4 + CH 3 




18.4 6 


16.6 C 


15. l a 


17.0 


16.1 


17.4 


15.5 


15.7 


13.1 


14.5 


38 


3 


CH 3 F + H -» CH 2 F + H 2 




14.2 ft 


12.9 C 


10. 4 a 


13.3 


14.6 


11.5 


10.2 


7.5 


10.9 


11.8 


39 


3 


NH 3 + H -> NH 2 + H 2 




16.3 ft 


14.9 C 


10. l a 


15.7 


15.9 


14.1 


13.2 


9.0 


12.7 


13.9 


40 


3 


NH 3 + O -> NH 2 + OH 




16.7 6 


12.2 C 


11. 4 a 


10.6 


9.8 


12.5 


11.8 


4.3 


4.0 


6.5 


41 


3 


CH 4 + 2 -» CH 3 + H0 2 




55.6 i > 


56.7 C 


56. C 


56.2 


55.0 


55.2 


54.4 


52.8 


53.3 


56.6 


42 


3 


CH 3 F + O -+ CH 2 F + OH 




14.3 fc 


11. 7 C 


11. a 


9.4 


9.8 


11.6 


10.3 


4.8 


3.9 


5.8 


43 


3 


CH 3 F + OH ->■ CH 2 F + H 2 




6.7 6 




4 2 a 


3.3 


5.0 


5.3 


3.2 


-0.0 


-1.1 


0.3 


44 


3 


CH 2 F 2 + O ->■ CHF 2 + OH 




15.2 6 




11. 4 a 


9.9 


10.4 


12.1 


11.0 


5.0 


4.0 


5.9 


45 


3 


CHF 3 + OH ->■ CF 3 + H 2 




7.5 b 




7.8 a 


6.7 


8.4 


8.6 


6.8 


2.9 


1.5 


2.8 


46 


3 


CH 3 Br + O ->■ CH 2 Br + OH 




14.4 6 




12. l a 


10.2 


10.4 


12.5 


11.2 


5.2 


4.4 


6.3 


47 


3 


CHF 3 + O -> CF 3 + OH 




23.9 ft 




15. 9 a 


14.6 


14.7 


16.6 


15.7 


9.4 


8.3 


10.2 


48 


3 


CH 2 F 2 + OH -» CHF 2 + H 2 




6.1 6 




4.8 a 


3.4 


5.1 


5.3 


3.4 


-0.1 


-1.2 


0.1 


49 


3 


CF 3 + H 2 ->• CHF 3 + H 




10.9 ft 




9.4 a 


11.0 


9.6 


10.3 


7.6 


9.4 


8.5 


9.9 


50 


3 


CH 3 Br + OH -» CH 2 Br + H 2 




6.2 6 




5.1 a 


4.0 


5.6 


6.0 


3.8 


0.5 


-0.6 


0.8 


51 


3 


CH 2 Br + H 2 -> CH 3 Br + H 




14.2 6 




11. 5 a 


13.1 


11.8 


12.6 


10.4 


11.9 


10.7 


12.3 


52 


3 


C 2 Hg + CH 3 — > C 2 H5 + CH 4 




16.0 ft 




12. 8 a 


15.1 


14.3 


15.4 


13.4 


13.7 


11.0 


12.3 


53 


3 


SiH 4 + H -> SiH 3 + H 2 




5.6 6 




3.0 a 


5.6 


6.5 


3.6 


1.9 


1.3 


5.1 


5.8 


54 


3 


NH 3 + NH 2 -> NH 2 + NH 3 




14.6 6 


13.7 C 




10.9 


10.9 


12.0 


10.5 


7.0 


5.2 


7.0 


55 


3 


PPH+ -> HP-PH+ 




39.3 ft 






41.9 


41.3 


41.2 


42.4 


39.2 


38.9 


37.9 


56 


3 


CH 4 + N0 2 -> CH 3 + HN0 2 




35.9 6 






36.9 


36.7 


34.5 


32.5 


34.3 


33.3 


35.1 


57 


3 


H 3 CO -► H 2 COH 




32.9 6 


31.4 C 




34.5 


36.8 


35.6 


36.2 


33.5 


32.9 


32.6 


58 


3 


H0 2 + NH -» 2 + NH 2 




1.9 6 


0.6 C 




-0.0 


0.2 


1.2 


0.1 


-4.4 


-5.5 


-4.0 


59 


3 


CH 2 F 2 + NH -> CHF 2 + NH 2 




21. 7 b 


13.0 C 




18.0 


17.2 


18.5 


17.2 


13.9 


12.2 


14.0 


60 


3 


NH 2 + O -> NH + OH 




9.9 b 


6.3 C 




4.7 


5.8 


7.3 


5.8 


-0.5 


-1.4 


0.1 


61 


3 


CH 3 + O -> CH 2 + OH 




n.o b 


14.9 C 




13.4 


15.6 


15.4 


14.0 


10.0 


8.6 


10.2 


62 


3 


CH 3 CH 2 CH 2 CH 3 




42.8 b 


42.2 C 




42.6 


42.3 


43.0 


43.1 


41.5 


40.7 


40.3 


63 


3 


CH 3 F + H -> CH 3 + HF 




31.4 6 


30. 9 C 


30. 6 a 


31.4 


31.1 


30.5 


30.4 


21.7 


27.4 


29.5 


64 


3 


CF 4 + H — > CF 3 + HF 




40.9 6 




44.3 a 


44.2 


44.2 


44.2 


46.0 


32.7 


38.5 


40.7 


65 


3 


CH 3 Br + H ->■ CH 3 + HBr 




6.1 b 




5.9 C 


6.5 


7.7 


5.1 


3.4 


2.0 


5.7 


7.1 


66 


3 


so + o 2 -» so 2 + O 




18.3 6 




4.6 a 


6.4 


5.0 


10.3 


5.7 


6.3 


5.0 


8.4 


67 


3 


CO + o 2 -► co 2 + O 




52.3 b 




47.2 a 


51.8 


51.7 


48.1 


47.4 


47.2 


47.1 


49.7 


68 


3 


SCO + O — SO + CO 




2.4 6 




5.4 a 


3.0 


2.9 


5.9 


3.8 


-1.6 


-2.5 


-0.4 


69 


3 


CO + BrO -+ C0 2 + Br 




11.4 ft 




7.4 a 


12.8 


13.3 


12.4 


10.5 


9.6 


8.6 


10.2 


70 


3 


N 2 + H — > N 2 H 




14.5 6 


14.1 c 




13.3 


14.6 


11.2 


9.5 


8.2 


12.0 


13.2 


71 


3 


N 2 + Br -> N 2 + BrO 








33.0 a 


31.8 


31.0 


35.3 


32.2 


30.4 


29.1 


34.2 



a Energy has been corrected by scaled (0.9804) B3LYP/6-31G* ZPVE to obtain bottom-of the well values. 

b Energy has been corrected by scaled (0.8929) HF/6-31G* ZPVE to obtain bottom-of the well values. 
c Energy has been corrected by scaled (0.92511) HF/6-311G** ZPVE to obtain bottom-of the well values. 

d At mPWlK/6-31+G* geometry. 
e Reactant at B3LYP/6-31+G* geometry, transition state at mPWlK/6-31+G* geometry. 
f Using a 6-311+G(2d,p) basis set to optimize the geometries. 
* Dissociated to the products. 



